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ABSTRACT 

We have determined a dynamical mass for the companion to HD 33636 which 
indicates it is a low-mass star instead of an exoplanet. Our result is based on 
an analysis of Hubble Space Telescope {HST) astrometry and ground-based ra- 
dial velocity data. We have obtained high-cadence radial velocity measurements 
spanning 1.3 years of HD 33636 with the Hobby- Eberly Telescope at McDonald 
Observatory. We combined these data with previously published velocities to 
create a data set that spans nine years. We used this data set to search for, and 
place mass limits on, the existence of additional companions in the HD 33636 
system. Our high-precision astrometric observations of the system with the HST 
Fine Guidance Sensor Ir span 1.2 years. We simultaneously modeled the radial 
velocity and astrometry data to determine the parallax, proper motion, and per- 
turbation orbit parameters of HD 33636. Our derived parallax, TTabs = 35.6 ± 
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0.2 mas, agrees within the uncertainties with the Hipparcos value. We find a 
perturbation period P = 2117.3 ± 0.8 days, semimajor axis = 14.2 ± 0.2 
mas, and system inchnation i = 4?1 ± 0?1. Assuming the mass of the primary 
star Ma = 1.02 ± 0.03 Mq, we obtain a companion mass Mb = 142 ± 11 Mj^p 
= 0.14 ± 0.01 Mq. The much larger true mass of the companion relative to 
its minimum mass estimated from the spectroscopic orbit parameters (M sin i = 
9.3 Mjup) is due to the near face-on orbit orientation. This result demonstrates 
the value of follow-up astrometric observations to determine the true masses of 
exoplanet candidates detected with the radial velocity method. 

Subject headings: astrometry - planetary systems - stars: individual (HD 33636) 
and distances - techniques: radial velocities and interferometric 



1. INTRODUCTION 

Mass is the most important physical property of exoplanets because it is the one pa- 
rameter that can establish a star's companion as an actual planet rather than a brown 
dwarf or low-mass star. Exoplanet masses also provide important constraints for theoretical 
models. As for stars, mass critically determines most of the instantaneous characteristics 
and long-term evolution of a planet. Therefore, we can test, and ultimately improve, our 
understanding of planet formation and evolution by comparing model predictions with the 
observed physical and orbital properties of exoplanets with measured masses. 

Currently, fewer than 10% of the more than 200 candidate exoplanet^ orbiting nearby 
stars have precisely determined masses. The rest only have known minimum masses. This is 
because the most successful technique for detecting candidate exoplanets, the radial velocity 
method, cannot be used to remove the degeneracy between the mass and orbital inclination 
for most of the known exoplanet candidates. In principle, radial velocities alone can be 
used to determine the masses of exoplanets in multi-planet systems where two or more 



Nauenbere 


2002; 


Rivera et al. 


2005) 



(e.g. 

876), results from different groups vary significantly, and the effects of non-coplanarity have 
yet to be considered. Thus, establishing precise masses, rather than arguing statistically, for 
the majority of exoplanet candidates requires observations with complementary techniques. 
The techniques that have been employed to break the mass - inclination degeneracy in 



A regularly updated list of reported exoplanets can be found at 



http://exoplanet.eu/ 
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radia l velocity data are astrometry (e.g. [Benedict et al.ll2002al ) and transit (e.g. iHenry et al. 
20001 ) observations. In this paper we report our determination of a dynamical mass for the 
exoplanet candidate companion to HD 33636 based on combined modeling of radial velocity 
and astro metric measurements of its perturbation orbit. 



H D 33636 (= G 97-2 5, HIP 24205) is GO V star at a distance of 28.7 pc flPerrvman et al. 



19971 ). IVogt et al.l (120021 ) reported the discovery of a candidate planetary mass companion 
orbiting HD 33636 in a long p eriod and ecc e ntric orbit based on its radial velocity variations. 
Perrier et allJ2003h . and later lSutler et all J2006h. refined the c ompanions's orbit parameters 
based on additional velocity measurements. iButler et al.l (120061 ) found the spectroscopic orbit 
parameters period P = 2128 days, eccentricity e = 0.48, and velocity semiamplitude K = 164 
^ which resulted in an estimated minimum mass M sin i = 9.3 Mj^p for the companion. 



m s 



Assuming the iButler et al.l (120061 ) orbit parameters, the corresponding minimum astro- 
metric perturbation, qa sini, of HD 33636 due to its low-mass companion would be 1.0 mas. 
This suggested that the perturbation from even an edge-on system orientation {i = 90°) 
would be detectable using the Hubble Space Telescope ( HST) F ine Guidance Sensors (F GS), 
and was the ni o tivati on for this work. [Benedict et al.l (j2002al ). iMcArthur et al.l (120041 ). and 
Benedict et al.l (120061 ) have previously used the HST FGS in combination with radial veloc- 
ity measurements obtained with ground-based telescopes to directly determine the masses 
of three exoplanets. Our observations and analysis to determine the mass of HD 33636 's 
companion are similar and described in the following sections. 

In §2 we describe our radial velocity measurements of HD 33636 with the Hobby- Eberly 
Telescope (HET), which we have used to supplement the previously published velocity data 
for this object and search for additional low-mass, short-period companions. In §3 we de- 
scribe our HST observations and astrometry measurements. Our reference star spectroscopy, 
photometry, resulting spectrophotometric parallaxes, and estimated proper motions are dis- 
cussed in §4. In §5 we present our simultaneous modeling of the radial velocity and astrom- 
etry data. From this modeling we determine the perturbation orbit parameters, which allow 
us to calculate the mass of the companion to HD 33636. In §6 we show our analysis of the 
radial velocity residuals from the single companion model and present a calculation of our 
detection limits for additional companions. We discuss the implications of our result in §7 
and conclude in §8 with a summary. 
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2. RADIAL VELOCITY DATA 

The ultimate precision and accuracy of our measurement of the mass for HD 33636 's 
unseen companion depends strongly on the quality and quantity of the radial velocity data 
used. Unaccounted for, but detectable, additional companions or poorly constrained spec- 
troscopic orbit parameters for the known companion would each cause a systematic error 
in our result. We therefore carried out high-cadence spectroscopic observations for radial 
velocity measurements of HD 33636 with the HET. We used these data to supplement the 
previously published velocities for this object and the combined data set of radial velocities 
spans just over nine years. 



2.1. HET Spectroscopic Observations 

We obse rved HP 33 636 on 65 nights using the HET to feed the High Resolution Spectro- 



graph (HRS,llul]Jll998|) between UT dates September 20, 2005 and January 21, 2007. The 
HET is queue, rather than classically, scheduled, and the nights were randomly distributed 
throughout the observing seasons that the star was available. The HRS was used in the res- 
olution R = 60,000 mode with a 316 gr mm~^ echelle grating. The cross- dispersion grating 
was positioned so that the central wavelength of the order that fell in the break between the 
two CCD chips was 5936 A. A temperature controlled (T = 70?0±0?1 C) iodine cell (HRS3) 
was inserted in front of the spectrograph slit entrance for all exposures to imprint lines that 
provided a contemporaneous wavelength scale and instrumental profile (IP) fiducial. Three 
separate exposures were taken within 15 minutes on all but a few nights. In total, 195 "tar- 
get" observations of HD 33636 were made, including 64 sets of three observations each and 
three additional solitary exposures. Exposure times were nominally 120 s, but varied up to 
twice that occasionally to account for increased seeing and/or cloud cover. 

Additionally, we observed HD 33636 once on December 6, 2006 without the iodine cell 
and with the same instrument setup, but in the R = 120,000 mode. The exposure time for 
that observation was 600 s. We used the spectrum from this "template" observation as the 
model template in the radial velocity analysis of the target spectra described in §2.2. 

CCD reduction and optimal ord er extraction were carried out for all the individual 



spectra using the REDUCE package (jPiskunov fc Valenti 1 120021 ). The final median signal 



to- noise (S/N) per pixel of the target spectra is 89 and the template spectrum S/N per pixel 
is 188. For the CCD binning we used (two pixels in the cross-dispersion direction and one 
in the echelle dispersion direction), there are roughly four pixels per resolution element in 
the R = 60,000 mode spectra and two in the R = 120,000 spectrum. 
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2.2. Radial Velocity Analysis 

We us ed an independent a dapt ation of the canoni cal spectrum modeling technique de- 



scribed by IValenti et al.l (119951 ) and iButler et al.l (Il996l ) to measure relative radial velocities 



in the extracted target spectra. The observed spectra between 5020 and 5860 A were broken 
up into 637 separate 100 pixel "chunks." Each chunk was modeled as the product of a high- 
resolution FTS spectrum of the iodine celj^ and a template stellar spectrum convolved with 
an IP. The template spectrum was Doppler shifted before being multiplied by the iodine 
spectrum, and this constituted the measurement of the star's relative radial velocity. Addi- 
tionally, the wavelength scale of the observations was simultaneously modeled as a second 
order polynomial. The reduced observed spectra have an arbitrary normalization and the 
model normalization was a free parameter as well. 

We used the sum of 11 Gaussians, one central and five satellites on each side, to represent 
the HRS IP. The variable parameters were the width, cr, of the central Gaussian and the 
heights of the satellites. The height of the central Gaussian was calculated from the a with 
the standard formula {h = l/27rcr). The IP of the HRS varies considerably in shape and 
size along each echelle order. Therefore, we used a dynamic spacing of the satellites from 
the center of the central Gaussian. The spacing was set so that there were four satellites 
per the FHWM of the central Gaussian (P ~ 2.355a). The a of the satellites were also 
set dynamically. They were equal to the a of the central Gaussian multiplied by 0.35. 
Including the Doppler shift, three parameter wavelength scale, 11 parameter IP description, 
and spectrum normalization there was a total of 16 free parameters in the model for each 
spectrum chunk. 

We based our model template spectrum on the R = 120k HD 33636 spectrum taken 
without the iodine cell. We modeled a fiat field spectrum with the iodine cell that was 
obtained using the same instrument setup as the template observation, taken immediately 
after it in the same manner as for the velocity measurements, but without the Doppler 
shift parameter (15 free parameters). This yielded an accurate wavelength scale and IP 
description as a function of the spectral position on the CCD at that time. We used this 
information to estimate the intrinsic spectrum of HD 33636 by removing the IP f rom the 



template observation. We used a modified Jansson technique (iGilliland et al.lll992l ) for the 
deconvolution. This estimate of the HD 33636 intrinsic spectrum was used as the model 
template for the velocity measurement of the target spectra. 



^The FTS spectrum of the HRS3 iodine cell is available at 
ftp://nsokp.nso.edu/FTS-cdrom/FTS50/ 001023R0.004l 
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The spectral modeling yielded a velocity measurement for each of the 637 chunks in a 
spectrum. For each target spectrum, we calculated the weighted mean of the velocities from 
the chunks in a spectrum to give a single measured velocity and uncertainty. We weighted 
the velocities for each chunk according to the velocity error computed acco rding to the 



formula for the intrinsic Doppler error in a section of stellar spectrum given by iButler et al. 



( 119961 ). The median uncertainty in our 195 measured weighted mean velocities is 2.1 m s~^. 
We converted the velocity determined from each spectrum into a relative radial velocity by 
correcting each measurement for the barycentric motion of the observatory in the line-of-sight 



direction using the JPL ephemeris DE40a^i. 

The RMS deviation of the HET velocities from the perturbation orbit we determine 



see §5) is 6.1 m s ^. Using the method of I Wright I (120051 ) . iButler et al.l (120061 ) estimated 



the intrinsic radial velocity "jitter" due to variations of the stellar photosphere for HD 33636 
to be 5.2 m s^^. Adding this value in quadrature to our velocity uncertainties we find xt 
= 1.1 for the fit to the orbit, which indicates general agreement between our radial velocity 
measurement and error estimation technique and that of the California - Carnegie Planet 
Search group. 

Our scheme of taking three successive exposures over the course of ~ 10 minutes during 
most nights gives us the opportunity to somewhat reduce the impact of short-term stellar 
noise and random errors on our orbit ana lysis. In addition to calculating a simple weighted 



mean, we also used the program Gaussfit ( iJefferys et al.lll988l ) to calculate a robust weighted 



midpoint time, velocity, and uncertainty for each of the 64 three observation sets. Comparing 
the resulting data to the fit for the weighted and robust mean methods, we find RMSs of 4.4 
and 3.3 m s~^ respectively. Because of the higher quality of the robust combined velocities, 
we ultimately adopted those data in the analysis presented in §5. The velocity values, 
velocity uncertainties, and the corresponding heliocentric Julian dates are given in Table [TJ 
The velocities we measured are relative to an arbitrary zero-point, and we have subtracted 
a constant value determined during our orbit analysis. Therefore, the velocities given are 
relative to the HD 33636 system barycenter under the assumption that our single-companion 
model for the system is correct. 



-"^^The JPL ephemeris data may be obtained on the Internet at 



ftp: / /ssd.jpl.nasa.gov / pub / eph/export /DE405 
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2.3. Total Radial Velocity Data Set 



In addition to our own measurements from the HET, we included pu blistied high preci- 
sion radi^J__v;e]^cityjii^suremen^ of HD 33636 from iPerrier et al.l (l2003l . "Elodie" sample) 
and iButler et al.l (l2006l . "Lick" and "Keck" samples). G. Marcy (pri vate communication ) 
provided us with an updated version of the Keck velocities published in iButler et al.l (120061 ) . 
These updated velocities result from a reanalysis of old spectra with an improved version of 
the California - Carnegie Planet Search Doppler analysis software and include an additional 
epoch of data taken since publication that overlaps with our HET observations. Table [2] 
lists the source, time coverage, number of data points, and RMS deviation from the final 
orbit fit for the individual samples. 

The complete data set contains radial velocities from four different telescope and in- 
strument combinations and has a time baseline of 3289 days (9.01 years). The preexisting 
data were a crucial component of our analysis, because HD 33636 's companion has an orbital 
period more than four times the time span of our HET observations. However, the use of 
a heterogeneous data set does require particular care, because the velocities in each sample 
are relative to a different zero point. We found in the course of our analysis (see §5) that a 
simple offset parameter for each sample was sufficient to correct the data to the same frame 
of reference. 



3. HST ASTROMETRY DATA 



We used the HST Fine Guidance Sensor Ir (FGSlr) to carry out astrometric observa- 
tions of HD 33636 and five reference stars between UT dates August 20, 2005 and October 



26, 20 06. A detailed overview of the FGSlr as a science instrument was given bv iNelan et al. 



(2003). Our data acquisition and reduction follow the procedure outlined by [Benedict et al. 



(I2OOOI ) as for the FGS3. We used the FGSlr for the current study beca use it provides superior 
fringes from which to obtain stellar positions (IMcArthur et al.ll2002l ). 



Table [3] lists the log of the HST observations. The data span 432 days (1.18 years) 
and include 18 epochs. Each epoch contains 2-4 positional measurements of HD 33636 
and the references stars, which were acquired contiguously over a time span of 26 - 35 
minutes. The observation time listed is for the midpoint of each epoch. The field was 
observed at multiple spacecraft roll values, and HD 33636 had to be placed in different non- 
central locations within the field of view (FOV) to accommodate the distribution of reference 
stars. The F5ND neutral density filter was used for the observations of HD 33636, while the 
F583W filter was used for the reference stars due to their being much fainter. To account 
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for using a different fil ter for the reference stars, we included a cross-filter correction term 
( iBenedict et al.ll2002bl ) in the astro metry model (see §5) for HD 33636. 



ASTROMETRIC REFERENCE STAR DATA 



Because of the high sensitivity of the FGS as an astrometer, and despite the relatively 
large distance of the astrometric reference stars, we had to take into account their parallaxes 
and proper motions in our model (see §5). To establish these, we determined spectropho- 
tometric parallaxes, and adopted proper motions that could be input as constraints to our 
model. Additionally, our model requires input [B — V) colors for all the stars to correct for 
chromatic aberration (lateral color). Our m ethod for estimating re ference star spectrophoto- 
metric parallaxes is discussed extensively in [Benedict et al.l (120071 ) and we followed the same 
approach in the present study. 

We obtained classification spectra of our astrometric reference stars with the R-C Spec- 
trograph on the Cerro Tololo Inter-American Observatory (CTIO) Blanco 4 m telescopj"]. 
The spectral types were determined from these spectra by a combination of template match- 
ing and line ratios and are generally better than ± 2 subclasses. 

We also obtained BVI photometry for the reference stars using the New Mexico State 
University Im telescope and adopt ed JHK magnitud es from the Two Micron All-Sky Survey 
(2MASS) Point Source Catalog (ICutri et al. I l2003l ). We estimated the star's luminosity 
classes with this photometry and the reduced proper motion method. We estimated V band 
extinctions. Ay, for the reference stars b y comparing their measured colors with the expected 
colors for their spectral type taken from I Cox I (2000). 



We next estimated My values for the reference stars by ass u ming the prototypical values 
for their spectral types and luminosity classes given by I Cox I (120001 ). We calculated their 
spectrophotometric parallaxes from their measured V and the estimated Ay and My values. 
The determined spectral types and luminosity classes, measured V and [B — V) values, esti- 
mated Ay values, assumed My values, and estimated spectrophotometric parallaxes for the 
reference stars are given in Table HI Proper motions for these stars were taken frora the 2nd 
data release of the USNO CCD Astrograph Catalog (UCAC2. IZacharias et al.ll2004l ). These 
proper motions, along with the spectrophotometric parallaxes and measured {B — V) values. 



^^CTIO is operated by the Association of Universities for Research in Astronomy Inc. (AURA), un- 
der a cooperative agreement with the National Science Foundation (NSF) as part of the National Optical 
Astronomy Observatories (NOAO). 
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enter our model as observations with error. 



SIMULTANEOUS RADIAL VELOCITY AND ASTROMETRY 

SOLUTION 



We modeled the radial velocity and astrometry data simultaneously to determine the 
parallax, proper motion, and complete set of perturbation orbit parameters for HD 33636. 
With these determined parameters, we then calculated the mass of its compani on. The 



method we used is ve ry si milar to that pr e viousl y employed by [Benedict et al.l (j2002al ). 
McArthur et al.l (12004 ) , and [Benedict et al.l (120061 ) to determine the same parameters for 
other exoplanet host stars and their companions. 

The astrometric model for HD 33636 that we used is represented by four solved equations 
of condition. They are 



x' = x + LC,{B-V)-XF,, (1) 

y' = y + LCy{B-V)-XFy, (2) 

i = Ax' + By' + C - P^n - fi^At - ORBIT^, (3) 

r] = -Bx' + Ay' + F- Pstx - fisAt - ORBITs. (4) 



Identifying terms, x and y are the measured coordinates; (B — V) is the photometric color; 
LCx and LCy are the lateral color corrections; XF^ and XFy are the cross-filter corrections; 
A and B are plate scale and rotation parameters; C and F are offsets; fia and fis are proper 
motion components. At is the time difference from the mean epoch; and Ps are parallax 
factor components; vr is the parallax; ORBITs and ORBITs are the astrometric components 
of the perturbation orbit; and ^ and r] are the standard coordinates. The astrometric orbit is 
a function of the orbital parameters period (P), time of periastron passage (Tp), eccentricity 
(e), semimajor axis {a a)-, position angle of the ascending node (f2), inclination (i), and the 
longitude of periastron passage {u). 

We also modeled the astrometry for the five reference stars in parallel with HD 33636. 
The equations of condition for those stars are the same as given in Equations (1-4) minus 
the cross-filter corrections and perturbation orbit motion. The plate scale and rotation 
parameters and the offsets are the same for HD 33636 and the reference stars at each epoch. 
The parallax and proper motion are unique for each star. 

The radial velocity model for HD 33636 is given by a fifth solved equation of condition, 

^ = RV + Gs- ORBITr, (5) 
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where RV is the measured relative radial velocity; Gs is the velocity offset for each of the 
four velocity samples described in §2.3; ORBITji is the radial component of the orbital 
velocity; and 7 is the adopted velocity of the HD 33636 system barycenter. The exact choice 
of 7 is arbitrary and immaterial because our analysis is based on the relative motion of HD 
33636. The important point is that its value is the same for each radial velocity sample and 
the Gs values are used to correct the sample velocities to the same frame of reference. The 
radial velocity orbit depends on the parameters P, Tp, e, and cj, which are the same for the 
astrometric and radial velocity orbit models. The radial velocity orbit is also dependent on 
the velocity semiamplitude {Ka)-, which does not influence the astrometric orbit. 

In addition to the shared orbit parameters, we enforce d a relationship betw e en th e 



astrometric and radial velocity models using a constraint from iPourbaix fc Jorissen I (120001 ). 



a A sin i _ PKa\/1 - 
'Kahs ~ 27r(4.7405) ' ^ ' 

where astrometric only quantities are on the left. Quantities determined primarily or only 
by the radial velocities are on the right. 

We solved for HD 33636's cross filter correction, coordinates, parallax, proper motion, 
and perturbation orbit parameters, the reference stars' coordinates, parallaxes and proper 
motions, the plate parameters for each astrometry epoch, lateral color corrections, and the 
velocity sample offsets by fitting the above models to the radial velocity and astrometry 
data. The HD 33636 cross-filter correction, lateral color corrections, and the reference star 
spectrophotometric parallaxes and proper motions were input into the model as observations 
with error. We used the Gaussfit program with robust estimation and the "fair" metric 
to determine the parameter values that gave the lowest between our model and the 
measured data. We adopted the uncertainties returned by Gaussfit, which were generated 
from a maximum likelihood estimatio n that is an ap proximation to a Bayesian maximum a 



posteriori estimator with a flat prior (jjefferys Ill99 



app 



The radial velocity data, our best fit, and the fit residuals as a function of time are 
plotted in Figure [H The radial velocities for each sample in the Figure have been corrected 
to the system barycenter using the corresponding offset values determined in the analysis. 
The RMS deviations for the individual samples are given in Table |2l We find RMS residuals 
for the Keck and HET velocities of 4.2 and 3.3 m s~^ respectively. 

The right ascension and declination component astrometry data as a function of time 
with the plate scale variation, parallactic motion, and proper motion removed are shown in 
Figure [2l These data illustrate the astrometric orbital motion of HD 33636. The fit to the 
data is indicated and the RMS residuals from our fit for HD 33636 are 1.4 mas and 1.0 mas in 
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the right ascension and dechnation directions respectively. These values are consistent with 
the previously characterized FGSlr per observation precision of ~ 1 mas and the median 
residual RMSs for the 5 reference stars, which is 1.3 mas. The quality of the simultaneous 
orbit fit to the astrometry and radial velocity data is evidence of the clear detection of HD 
33636's perturbation due to the known companion. The same data from Figure [2] are also 
shown as in the fiat plane of the sky in Figure [31 The direction of orbital motion and the 
location and time of the next periastron passage (2010.65) are also indicated in Figure [31 

Relative coordinates and our determined absolute parallaxes and proper motions for HD 
33636 and the five reference stars are given in Table [51 A summary of the HST astrometry 
for HD 33636 is given in Table [61 For HD 33636 we find Hats = 35.6 ± 0.2 mas, /Iq, = 
169.0 ± 0.3 mas, and fis = -142.3 ± 0.3 mas. Our parallax value is in exc ellent agreement 



with, but more precise than, the Hipparcos value iiabs = 34.9 ± 1.3 mas (jPerryman et al. 



19971 ). However, our determined proper motion is signifi cantly different than t he Hipparcos 



values /ia = 180.8 ±1.1 mas and fis = -137.3 ± 0.8 mas (jPerryman et al.lll997l ). The reason 



for this is that the Hipp arcos model did n ot ac count for the large perturbation orbit of HD 



33636. As discussed by [Black &: Scargld (1l982l ). the proper motion parameters can absorb 



orbital motion in astrometry data. The Hipparcos satellite only made 16 measurements of 
HD 33636's position with a reported precision of better than 5 mas over a period of 2.1 years. 
The orbital motion of HD 33636 during this time led to the underestimation of southerly 
proper motion component and overestimation of the westerly proper motion component. 
The Hipparcos measurements were not precise enough and did not span a sufficient amount 
of time to distinguish between proper motion and orbit curvat ure at the level necessary to 



trigger a multiphcity flag in the standard analysis (lESA |[l997l ) 



In the case of our analysis, we have the benefit of foreknowledge about the existence 
of a companion due to the radial velocity variations. The radial velocity data that we 
modeled simultaneously with the astrometry data is very sensitive to most of the orbit 
parameters and carries the most weight in the determination of those parameters because 
of its quantity and quality. Additionally, we have five times the per observation astrometric 
precision and four times the number of observations than Hipparcos for HD 33636. The 
Hipparcos data are not useful to include in our analysis because that would necessitate the 
inclusion of additional transformation terms similar to the velocity offsets used to combine 
the radial velocity datasets. This would result in further degradation of the Hipparcos 
precisions because there is no astrometric plate overlap between our HST observations and 
those of Hipparcos. Because of this, it would not have been useful to include the Hipparcos 
measurements in our analysis. 

Our derived values and uncertainties for the perturbation orbit parameters are given in 
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Table [71 Figure JH shows a plot of the r elatio nship between aA_ and i for fixed P, Ka, e, and 
TTabs through the iPourbaix fc Jorissen I (120001 ) relationship (eq. [6]). We find that the orbit 
is nearly face-on, with z = 4?1 it 0?1 from the plane of the sky. Correspondingly, we find a 
large perturbation size, = 14.2 ± 0.2 mas, relative to the previously calculated minimum 
perturbation a^sini = 1.0 mas. We also find P = 21 17.3 ± 0.8 days. Assuming the mass 
of HD 33636 Ma = 1-02 ± O.O3M0 ( iTakeda et al.l 120071 ). we find the mass of the companion, 
Mb, by iterating the equation 



T^abs 



(7) 



This yields Mb = 142l^;8 Mjup = 0.136To;oo2 Mq, which is also relatively large compared to 
the minimum mass Msini = 9.3Mjup, calculated from the spectroscopic orbit parameters. 
We have elected to round up to the next significant figure in solar units to account for possible 
systematic errors in the adopted mass of HD 33636 and other aspects of the analysis. Our 
final adopted companion mass uncertainty, which is also given in Table [TJ is O.OIM© = 
llMjup- We conclude from this result that HD 33636's companion is a low-mass star and 
not an exoplanet. 



6. LIMITS ON ADDITIONAL COMPANIONS 



As discussed in §2, one motivation for obtaining additional radial velocity data beyond 
the previously published data was to search for, and place limits on, additional companions 
in the system. To do this we analyzed the radial velo city residuals from the orbit model 
described in §5. We calculated a (Lomb) periodogram (jPress et al.lll992l ) for the residuals 
from the model fit and the result is shown in Figure O We searched for prominent peaks 
at periods shorter than the time span of the data set (3289 days) and found none with 
significant power or having a < 70% false alarm probability (FAP). We also carried out this 
analysis on the residuals for the individual samples separately and, again, found no signs for 
regular periodicity. In addition, we analyzed the radial velocity data by fitting trends to the 
total data set and individual samples' residuals and found nothing significant. For the total 
data set we found a trend of 0.3 m s~^ over the 9 year span, which is within the velocity 
amplitude uncertainty (0.6 m s~^, see Table [7]). 

To q uantify our detection limit s for additional companions we used a method similar 
to that of IWittenmyer et al. I (120061 ). We adopted the velocity residuals as a noise sample 
for simulated orbits. We generated simulated radial velocity orbits for each of 200 period 
values from 0.8 to 3289 days. The discrete period values were spaced evenly on a logarithmic 
scale. For each period, we began by assuming a starting velocity semimajor amplitude of 
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4 m s~^. We calculated 100 orbits for the given period and velocity semiamplitude and 
selected different time and longitude of periastron passage values with a pseudo-random 
number generator. We calculated the radial velocity values of the orbits at each of the dates 
in our measured velocity data set. We scrambled our velocity residuals, also with a pseudo- 
random number generator, and added them to the simulated data as noise. For each of the 
100 trials we calculated a periodogram of the simulated velocity data. If the power at the 
period used to generate the orbit was found to correspond to a FAP < 0.1% then that was 
counted as a successful detection. If 99 of the 100 trials at a velocity semimajor amplitude 
resulted in detections then that velocity was considered the detection limit at the period 
value and the next period value was considered. If two or more non-detections occurred 
then the velocity semimajor amplitude was increased 0.5 m s~^ and the process continued 
until the velocity limit for the simulated period was found. We carried out this analysis for 
three different values of orbital eccentricity e = 0.0, 0.4, and 0.7. The upper eccentricity 
limit of 0.7 was chosen because > 95% of radial velocity detected exoplanet candidates have 
eccentricities lower than that value. For the simulations with an eccentricity of 0.0, the 
longitude of periastron is a meaningless parameter, and the 100 randomized values of the 
time of periastron served to set the relative phase. 

With the set of eccentricities, periods, and velocity detection limits derived from this 
procedure we calculated the corresponding minimum companion masses and astrometric 
perturbation sizes. The results are shown in Figure El From these data we infer that 
our radial velocity data were at least sensitive enough to detect companions with minimum 
masses Msini > 0.2 Mjup, orbital periods P < 100 days, and orbital eccentricities < 0.4. 
For the same orbital eccentricities, companions with Msini > 1.2 Mjup and orbital periods 
up to 3300 days are ruled out. We also find that radial velocity undetected companions 
would have a minimum astrometric signature asini < 0.05 mas for P < 1000 days and 
reasonable eccentricities, which is beyond the detection threshold of the HST observations. 
Taken together, the periodogram analysis, slope analysis, and simulation of the residuals 
indicate no additional companions in the system, verifying that our single companion model 
to fit the radial velocity and astrometry data was appropriate. 



DISCUSSION 



HD 33636's compan ion is likely a M6 V star, with magnitudes AV ~ 8.0 and AK ^ 4.0 
from the GO V primary (jCox Il2000l ). From this we estimate that HD 33636 should have only 
a 0.03 mag enhancement in K from that expected for a solitary GO V star. The two stars' 
separation at apastron will be 0'.'15. Therefore, high-resolution imaging could be useful for 
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additional study of this system. 



The mass we find for HD 33636's companion indicates that it is not a planet at the 
~ 12(7 level assuming the standard 13 Mj^p upper limit for planets. This is the first definitive 
example of a planet candidate detected with the ra dial velocity method beinR con firmed to 
exist, but to have a non-planetary mass. Previously. iReffert fc QuirrenbachI (120061 ) analyzed 
the Hipparcos Intermediate Astrometric Data for HD 38529 and 168443 with radial velocity 
data constraint and found companion masses M = 37tlt Mjup and M = 34 ± 12 Mj^p 
respectively. However, these are both < 1.8 cr results and more observations are needed to 
refine these estimates. 



In contrast, Benedict et all fl2002ar ). iMcArthur et al.l (120041 ) . and Benedict et al.l (120061 ) 
have directly confirmed the planetary nature of exoplanets around GJ 876, Cancri, and e 
Eridani with the same method used in this paper. The GJ 876 and Cancri systems contain 
additional, non-astrometrically detected companions. If these systems are coplanar then the 
additional companions are also planets. Three other planet candidates originally detected 
with the radial vel o city method have been ob served to transi t their host star, H D 209458b 
jHenrv et allboool: Icharbonneau et"aDl200oh . HD 189733b JSouchv et al.ll2005h . and HD 



149026b JSato et al.lboosl ). and thus have measured masses and are confirmed planets. 



From mathematical arguments, the median value for the inclination of binary (star + 
star or star + planet) orbits is 60°. This is supported by the distrib ution of the inclinations o f 
visual binary star orbits in the Washington Double Star Catalog (IMason &: Hartkopill2006l ). 
which has a broad peak in frequency centered around 60°. An inclination i = 60° corresponds 
to a multiplicative factor of 1.15 to the minimum mass calculated from spectroscopic orbit 
parameters. Therefore, we expect that most of the candidate exoplanets detected with the 
radial velocity method and having minimum masses Msini < llMj^p are actually planets. 

Nearly face-on orbits, like the one we have determined for the HD 33636 system, should 
be rare. Orbits with i < 5° are expected to make up only 0.4% of an unbiased distribution. 
Nevertheless, our result is a striking example proving that minimum masses are not true 
masses and that not all of the planet candidates are actual planets. This demonstrates the 
value of follow-up astrometric observations and photometric monitoring for potential tran- 
sits to determine the true masses of exoplanet candidates detected with the radial velocity 
method. 



Valenti fc Fischer I (120051 ) found HD 3363 6's iron abundance [Fe/H] = -0.13, which is 



consistent with the solar neighborhood average (lAUende Prieto et al.ll2004j), but on the lower 
end o f the distribution of stars with detected planet candidate companions (IFischer fc Valenti 
20051 ) . Although HD 33636's companion falls outside the period range (P < 4 years) consid- 
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ered by iFischer fc Valenti I (120051 ) . it is reasonable to assume that the metalhcities of host 
stars to high-mass planets at all periods should be distributed at similar high values as is 
sug gested by the core accretio n model of planet formation. For stars with -0.50 < [Fe/H] < 
0.0, iFischer &: Valenti I (120051 ) found a planet candidate detection rate < 3%. Above solar 
metallicity, they found the detection rate increased, and was up to 25% for stars with [Fe/H] 
> -t-0.3 dex. In this context, the fact that HD 33636 is not a planet hosting star is not as 
surprising as it would have been if it had [Fe/H] ^ 0.0, because removing it from the sample 
strengthens the correlation between planets and host star metallicity. 



Along the lines of planet host star abundances, lEcuvillon et al.l (120061 ) and lChen &: Zhao 



(120061 ) have included HP 33636 as a planet ho sting star for studies of oxygen and lithium 
abundances respectively. lEcuvillon et al.l (120061 ) found that planet hosting stars could have 
oxygen abundances, [0/H], enhanced from a volume limited control sample by 0.1 - 0.2 dex, 
but that there was a large unc ertainty in the n ieasur ements and ambiguity with the effects 
of galactic chemical evolution. IChen &: Zhao I (120061 ) found that planet hosting stars could 
have depleted lithium abundances relative to stars without detected planet companions. HD 
33636 was found to b e slightly oxygen rich ( ^ 0.1 dex) relative to control samp l es sta rs with 
similar [Fe/H] in the lEcuvillon et al.l J2006h study. Conversely, IChen fc Zhao I (120061 ) found 
that HD 33636 showed no signs of lithium depletion as was the case in many other ostensibly 
planet hosting stars. 

Another st u dy w ith conclusions that could be affected by our result was done by 
Beichman et al.l (120051 ) . who looked for infrared excess due to debris disks around planet 
hosting stars with Spitzer. HD 33636 was one of six planet hosting stars that showed excess 
emission at 70 /z m and it also was one of the three that showed the most significant ex- 
cess. In addition, iBeichman et al.l (120051 ) found six stars that did not have detected planets 
showing the same type of excess. Moving HD 33636 from the planet hosting sample to the 
non-planet hosting sample has the effect of increasing the offset in the frequency distribution 
for 70 /xm excess between the two samples. Surprisingly, the non-planet hosting stars then 
have a higher infrared excess detection rate (7/60) than the planet hosting stars (6/80). This 
indicates there is no special correlation between planets and debris disks within the current 
planet and disk detection limits. Both this study, and the two mentioned above, depend 
on data that is difficult to measure. Resolving a true planet population from the sample of 
planet candidates with mass measurements will ultimately increase the impact of these and 
other similar studies and also permit new studies to be undertaken. 
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8. CONCLUSIONS 

We have determined HD 33636 's parallax, proper motion, and perturbation orbit param- 
eters from simultaneously modeling radial velocity and astrometry data. The total radial 
velocity data set spans 9.0 years and includes 1.3 years of data obtained by us with the 
Hobby- Eberly Telescope. The astrometry data spans 1.2 years and was obtained with the 
FGSlr instrument on the Hubble Space Telescope. 

Our derived parallax agrees with the Hipparcos value within the measurement uncer- 
tainties. We find a system inchnation i = 4?1 ±0?1, which indicates a near face-on orbit and 
a companion mass much larger than the minimum mass calculated from the spectroscopic 
orbit parameters. With our determined perturbation period P — 2117.3 ± 0.8 days and 
semimajor axis qa — 14.2 ± 0.2 mas and assumed mass for HD 33636 Ma — l.O2±O.O3M0, 
we obtain the mass of the companion Mb — 0.14 ± O.OIM0. We conclude that HD 33636's 
companion is a M6 V star rather than an exoplanet. 

We are grateful to G. Marcy for providing the unpublished Keck radial velocities. In 
addition, we acknowledge R. P. Butler, D. A. Fischer, J. Johnson, K. Peek, S. S. Vogt, 
and J. T. Wright for their contributions to measuring those velocities. We also thank the 
anonymous referee for a careful reading of the manuscript and constructive comments which 
improved this paper. JLB, BEM, and GFB acknowledge support from NASA GO-09407, GO- 
09969, GO-09971, GO-10103, GO-10610, and GO-10989 from the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555; and from JPL 1227563 MASSIF Key Project, 
Todd Henry, P.I.), administered by the Jet Propulsion Laboratory. 
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Table 1. HEX Radial Velocities for HD 33636 



HJD - 2450000.0 RV (m s"!) 



3633, 


,9377 


85.6 


± 3.2 


3646, 


,9084 


64.7 


± 3.3 


3653, 


,9013 


74.5 


± 3.5 


3663, 


,8716 


65.8 


± 2.9 


3666, 


,8409 


66.9 


± 3.2 


3668, 


,8329 


62.5 


± 3.4 


3676, 


,8395 


61.8 


± 3.5 


3678, 


,8104 


58.0 


± 2.9 


3680, 


,8050 


56.1 


± 3.0 


3682, 


,7969 


55.4 


± 3.2 


3683, 


,8119 


48.8 


± 3.3 


3689, 


,9219 


44.3 


± 3.2 


3691, 


,7897 


51.9 


± 3.0 


3692, 


,7895 


47.8 


± 2.9 


3696, 


,7711 


47.4 


± 3.0 


3697, 


,7683 


50.3 


± 2.9 


3700, 


,7600 


47.3 


± 2.6 


3703, 


,7523 


43.0 


± 4.0 


3708, 


,8624 


43.2 


± 3.3 


3709, 


,8785 


40.5 


± 3.4 


3713, 


,7238 


48.2 


± 3.9 


3714, 


,8699 


43.7 


± 3.9 


3719, 


,6983 


41.3 


± 4.1 


3719, 


,8439 


43.1 


± 3.9 


3724, 


,8191 


41.2 


± 3.4 


3724, 


,8225 


39.7 


± 3.6 


3730, 


,6675 


31.5 


± 3.7 


3731 


,6737 


40.8 


± 3.5 


3732, 


,6648 


41.4 


± 3.5 


3738, 


,6609 


38.7 


± 3.1 


3739, 


,6456 


32.5 


± 3.4 


3746, 


,6229 


37.7 


± 4.7 


3748, 


,6334 


30.4 


± 3.7 


3751, 


,7516 


25.9 


± 3.4 


3753, 


,7481 


28.9 


± 3.7 


3754, 


,6125 


28.6 


± 3.7 


3755, 


,6016 


28.7 


± 3.6 


3757, 


,7514 


24.6 


± 3.9 


3762, 


,5922 


26.2 


± 4.4 


3985, 


,9800 


-21.3 ± 3.4 


3987, 


,9655 


-31.4 ± 2.7 


3988, 


,9695 


-30.3 ± 2.7 


3989, 


,9691 


-31.3 ± 2.7 


3990, 


,9631 


-32.3 ± 2.7 


3997, 


,9516 


-27.1 ± 2.7 


4007, 


,9220 


-33.6 ± 2.8 
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Table 1 — Continued 



HJD - 2450000.0 RV (m s"!) 



4008.9051 


-32.9 


± 


3.1 


4014.9009 


-35.0 


± 


2.7 


4015.9059 


-33.7 


± 


3.3 


4018.8874 


-35.4 


± 


3.0 


4019.8780 


-36.6 


± 


2.8 


4020.8750 


-38.4 


± 


3.0 


4021.8729 


-37.5 


± 


3.0 


4031.8466 


-38.2 


± 


3.1 


4072.7382 


-45.0 


± 


3.4 


4073.7364 


-43.6 


± 


3.0 


4075.8628 


-46.1 


± 


3.0 


4076.7277 


-44.4 


± 


3.1 


4079.7194 


-37.6 


± 


3.3 


4080.8438 


-45.4 


± 


3.0 


4081.8592 


-48.9 


± 


3.2 


4105.6559 


-48.5 


± 


4.0 


4106.7734 


-50.6 


± 


3.7 


4108.7813 


-50.5 


± 


3.9 


4109.7746 


-51.8 


± 


3.7 


4110.7867 


-48.8 


± 


4.3 


4121.6098 


-51.6 


± 


4.1 
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Table 2. The Radial Velocity Samples 



Sample Time Span N RMS (m s"!) 



Lick 1998.05 - 2001.69 12 13.6 

Keck 1998.07 - 2006.68 27 4.2 

Ekxiic 1998.13 200.3.2.3 42 12.2 

llET 2003.72 2007.05 (i7 3.3 
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Table 3. Log of FGSlr Observations 



HJD - 2450000.0 




Hsr Roll (°) 


3603.2802 


4 


281.4 


3605.2786 


4 


281.0 


3606.5443 


4 


281.0 


3610.2767 


4 


281.0 


3613.1421 


3 


281.5 


3615.4771 


4 


281.0 


3658.0422 


2 


246.9 


3662.0402 


4 


247.6 


3667.3715 


4 


247.7 


3669.1033 


4 


250.9 


3670.2403 


4 


245.6 


3723.0689 


4 


182.0 


3967.5470 


4 


281.7 


3969.6785 


3 


281.0 


3972.3435 


3 


281.0 


4031.4272 


4 


251.9 


4033.1582 


4 


251.9 


4035.2237 


4 


251.9 



''Number of observations of HD 33636 
per epoch. 
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Table 4. Astrometric Reference Star Data 



Identification 


Spectral Type 


V 


B-V 


Av 




■Kabs (mas) 


Ref-1 


F6 V 


15.2 


0.6 


0.5 


3.6 


0.6 ± 0.1 


Ref-2 


F6 V 


14.1 


0.6 


0.4 


3.6 


0.9 ± 0.2 


Ref-3 


K6 V 


15.3 


1.3 


0.1 


7.3 


2.6 ± 0.5 


Ref-4 


G2 V 


13.1 


0.6 


0.1 


4.6 


2.0 ± 0.4 


Ref-5 


K3 III 


9.9 


1.3 


0.1 


0.3 


1.2 ± 0.2 



^Taken from lCo~ ] 1 I2OOOI) for the measured spectral types. 
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Table 5. Astrometry Catalog 



Star 




5=- 




Ma 






(arcsec) 


(aJTcsec) 


(mas) 


(mas) 


(mas) 


HD 33636 


-97.8347 ± 0.0003 


803.0143 ± 0.0005 


35.6 ± 0.2 


169.0 ± 0.3 


-142.3 ± 0.3 


Ref-1 


-172.1200 ± 0.0004 


800.5260 ± 0.0005 


5.7 ± 0.1 


8.6 ± 0.6 


-9.2 ± 0.5 


Ref-2 


-159.6280 ± 0.0002 


775.5080 ± 0.0005 


9.2 ± 0.1 


0.6 ± 0.3 


-12.2 ± 0.3 


Rcf-3 


-113.6041 ± 0.0004 


588.7119 ± 0.0005 


2.6 ± 0.2 


-4.1 ± 0.6 


4.8 ± 0.5 


Ref-4 


-216.5286 ± 0.0003 


618.6621 ± 0.0005 


2.0 ± 0.2 


-1.2 ± 0.4 


-2.5 ± 0.4 


Ref-5 


-262.7320 d= 0.0002 


583.6012 ± 0.0005 


1.2 ± 0.1 


0.0 ± 0.2 


-6.0 ± 0.2 


''The right 


ascension and declination coordinates are 


relative to a 


= 4''10'"22?66,<5 


= 5°13'29'.'0, 



J2000.0 
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Table 6. Summary of HST Astrometry 



Parameter Value 



HST study duration 


1.2 yr 


Number of HD 33636 obsrvations 


67 


Number of epochs 


18 


Number of reference stars 


5 


HD 33636 (V) 


7.1 


Reference stars {{V)) 


13.5 


HD 33636 (B-V) 


0.6 


Reference stars {{B — V)) 


0.9 


HST parallax* 


35.6 ± 0.2 mas 


Hipparcos parallax 


34.9 ± 1.3 mas 


HST proper motion^ 


220.9 ± 0.4 mas yr 


In position angle^ 


130? 1 d= 0?1 


Hipparcos proper motion 


227.0 ±1.4 mas yr 


In position angle 


127?2 ± 0?3 



''Prom the simultaneous modeling of the radial velocity and 
astrometry data. 
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Table 7. HD 33636 Perturbation Orbit Parameters and Companion Mass 



Parameter 


Value 


Ka 


163.5 ± 0.6 m s'^ 


P 


2117.3 ± 0.8 days 


To(JD) 


2451198.3 ± 2.0 


e 


0.48 ± 0.02 


U) 


337? ± 1?6 


a-A 


14.2 ± 0.2 mas 


n 


125? 6 ± 1?6 


i 


4?0 ± 0?1 


Mb 


142 ± 11 Mj^p"" 


Mb 


0.14 ± 0.01 Mq'^ 



'^Assuming = 1.02 ± O.O3M0 
jTakeda et al.ll2007l) . 



- 28 - 



200 7 




-20 - 



-40 



1998 2000 2002 2004 2006 2008 
Date [year] 



Fig. 1. — Upper panel: Radial velocities (points) as function of time and the determined 
orbit (dashed line) from the simultaneous fit to the radial velocity and astrometry data. The 
error bars are omitted for clarity. Bottom Panel: Residuals from the fit (points). 
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Fig. 2. — Right ascension (a) and declination (5) components of HD 33636's perturbation 
orbit as a function of time (points) and the best fit (dashed hne) from the simultaneous 
modeling of the radial velocities and astrometry. The dots are the individual observations 
and the circles are per epoch (single HST orbit) normal points consisting of 2 - 4 individual 
observations. 
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Fig. 3. — Perturbation orbit of HD 33636 on the sky (line). The open circles are the HST 
epoch normal points and are connected to the derived orbit by residual vectors. The HST 
data cover 20% of the orbit period. The orbital motion direction is indicated by the arrow. 
The square marks the location of periastron passage and the next time of occurrence is 
labeled. 
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Fig. 4. — The relationship (dotted hne) between t he perturbation size (a a) an d inchnation 
angle (i) for fixed P, Ka, e, and Habs through the IPourbaix fc Jorissen I (120001 ) relationship 
(eq. [6]). Our determined value for the perturbation size and inclination is given by the filled 
circle. The right axis maps the inclination to the corresponding companion mass (Mb)- Our 
adopted value for the uncertainty in the companion mass is plotted as the error bar for this 
axis. The formal uncertainties in our determined a a and i are smaller than the point. 
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Fig. 5. — Periodogram of the radial velocity residuals from the best fit, single-companion 
orbit with 1% and 10% false alarm probability (FAP) limits indicated. No periodicity is 
detected with FAP < 70%, which indicates that there are no additional companions in the 
system for our detection limits (see §6). 
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Fig. 6. — Results from the detection limit simulations converted to hypothetical companion 
minimum masses (M sin i, upper panel) and minimum astrometric perturbation size of HD 
33636 (oyisinz, lower panel). The different lines represent the different assumed eccentricity 
values. Values above the lines would have been detected with a periodogram analysis of the 
radial velocity data. 



